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■ Abstract 

A dynamical coupled-channel formalism for processes ttN — » KY and 7iV — » KY is presented 
which provides a comprehensive investigation of recent data on the jp — ► K + A reaction. The non- 
resonant interactions within the subspace KY (BttN are derived from effective Lagrangians, using a 
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unitary transformation method. The calculations of photoproduction amplitudes are simplified by 



casting the coupled-channel equations into a form such that the empirical jN — ► irN amplitudes 
are input and only the parameters associated with the KY channel are determined by performing 
X 2 -fits to all of the available data for ir~p — ► K°A, K°Y,° and ^yp — > K + A. Good agreement 
O ' between our models and those data are obtained. In the fits to irN — ► KY channels, most 

of the parameters are constrained within ±20% of the values given by the Particle Data Group 



X 



and/or quark model predictions, while for 7p — > i^ + A parameters, ranges compatible with broken 
SU(6) (g>0(3) symmetry are imposed. The main reaction mechanisms in i^ + A photoproduction are 
singled out and issues related to newly suggested resonances Sn, P13, and D13 are studied. Results 
illustrating the importance of using a coupled-channel treatment are reported. Meson cloud effects 
on the 7./V — > ^* transitions are also discussed. 
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I. INTRODUCTION 



lgh precision 
TO 

sjlaj covering 



Recent experiments at JLab-CLAS Q,Q, ELSA-SAPHIR 0, Q and Spring-8-LEPS Q, 
are refining our knowledge of associated strangeness photoproduction. 
differential cross section data for the process 7p — > fT + A have been released ^ 
the region between ll 7 pa 1.6 GeV and 2.6 GeV in the center-of-mass frame. Furthermore, 
Stag ,e po. arizatl0n as_ y data for recoil hype™ 2] and beam QQ have a, 8 o b e C o,„e 
available. 

The K + A photoproduction has also been extensively studied using phenomenological 



approaches. In general, those works 
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12j investigated the direct channel 
mechanisms based on an isobar approach in tree approximation. Combinations of isobar 



14. 
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models with a Regge analysis [13], successful at higher energies, have also focused 
on strangeness electromagnetic production. A new generation of more precise data has 
made it clear that coupled-channel effects can no longer be ignored and that multi-step 
processes have to be incorporated carefully. Coupled-channel formalisms based on the K- 



matrix approximation and isobar effective Lagrangians have been developed 



III 



The purpose of this work is to report on an advanced version of a dynamical coupled- 
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that incorporates proper treatment of off-shell effects. The 



channel formalism [13, 

direct KY photoproduction channel is investigated via a chiral constituent quark model 
(CQM) I22I 23 1 . This latter approach allows one to handle all known resonances with 
a reasonable number of adjustable parameters, in contrast to isobar effective Lagrangian 



models 24]. Consequently, the CQM provides an appropriate tool for understanding the 
elementary reaction mechanism, establishing reliable indicia for the predicted missing baryon 

30. 31. 



resonances [25, 



resonances. 
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321 ]. and gaining improved insights into the known 



In principle, the KY photoproduction should be investigated within a large scale coupled- 
channel approach including several reaction channels, e.g. ttN, rjN, uN, KY, <j)N, ttttN 
(criV, it A, pN). Obviously, this can not be done so easily because the data sets, to be 
simultaneously fitted, are very extensive, and reaction mechanisms involving channels other 
than 71N have not been studied extensively. 

As a first significant step, it is useful to consider a much more restricted coupled-channel 
model focusing on understanding particular reaction mechanisms. Concerning the KY pho- 
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toproduction, the obvious first task is to investigate the coupling between the KY and ttN 
channels for the following reasons. From the available data, one observes that kaon pho- 
toproduction is in general much weaker than pion photoproduction. Hence the multi-step 
transitions, such as 7iV — ► ttN — > KY, should be comparable to the direct 'yN — > KY 
process. This has been verified in Ref. using a coupled-channel model with 7 AT, ttN and 
KY channels. Moreover, the need for a coupled-channel approach to study meson-baryon 
reactions in the second and third N* regions has been well discussed in the literatures, as 



33 



34 ]. 



reviewed in Refs. 

In this work, we take advantage of the development of new models [19] for ttN — > KY 
and KY — > KY interactions to reinvestigate the influence of the ttN channel. Furthermore, 
we refine the models developed in Refs. jl^i Q] and consider recent 7p — > K + A data. 
Focusing on the coupled-channel effects associated with the ttN channel, we also determine 
the parameters of relevant N* resonances. Our results could serve as the starting point for 
performing more advanced coupled-channel calculations including additional meson-baryon 
channels. 

Within the considered coupled-channel model, a comprehensive study of K + A photopro- 
duction requires models of the non-resonant transitions among jN, ttN, and KY states and 
the decays into these three channels for about 12 isospin 1=1/2 N* states. In this work, 
wc follow Refs. HQ to derive the non-resonant transitions from effective Lagrangians by 
using a unitary transformation method [3] and SU(3) symmetry. For N* decays, we con- 
sider information from the Particle Data Group J^fl (PDG) and/or from constituent quark 
model predictions 22, Q- S] With these constraints, the model has a reasonable number 
of adjustable parameters, which can only be ascertained from the data. We simplify the 
fitting task by casting the coupled-channel equations into a form such that the empirical 
7iV — > TiN amplitudes [37[ are input to the calculations and only the parameters associ- 
ated with the KY channel are to be determined by performing \ 2 minimization fits to all 
available irN -> KY and 7N -> KY data, using the CERN-MINUIT code. 

In addition, to clarify the role of coupled-channel effects due to the irN channel, we 
also analyze the dynamical content of the jN — > N* transition. The so-called "meson 



cloud effects" discussed in the study 



of the A (1232) resonance are identified within 



our coupled-channel model. We also make an attempt to determine the properties o: 



predicted 



0, 22, 3, 0, 3 H H Q 



and/or sought for 
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14 L16, 20, 2J, |4 



41 



the 



42, 
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43|, |4J, |45|, |4y, |47[ third Su, Pis, and -D13 resonances. 

For simplicity, at this stage we do not consider the KT, photoproduction data to avoid 
the need to also determine the parameters associated with the photo-excitations of about 
12 other isospin 1 = 3/2 N* states. Obviously, our results serve as a good starting point for 
a subsequent investigation including all KY channels. Our results in that direction will be 
published elsewhere. 

This paper is organized as follows: In Section II, the theoretical frame is presented. The 
main content of our coupled-channel formalism is then given, followed by an outline of the 
relevant constituent quark model for the direct 7^ — > K + A channel. There, the novelties 
of our approach are discussed. Section III is devoted to numerical results and comparisons 
with available data for n^p — > K°A, K°Y,° and 7p — > K + A. For this latter reaction, the 
most relevant known nucleon resonances are singled out and possible manifestations of new 
baryon resonances are discussed. Meson cloud effects are exhibited by examining multipoles 
from the obtained model. Summary and conclusions are reported in Section IV. 

II. THEORETICAL FORMULATION 

In this Section, we first present our dynamical coupled- channel approach for the pho- 
toproduction process including intermediate ttN and KY channels. Then, we outline the 
constituent quark model used for the direct KY photoproduction reaction. 

A. Coupled-channel formalism 



The coupled- channe 



approach presented here is derived from a general formulation re- 



ported in Refs. |33l . |. The starting point is a Hamiltonian consisting of non-resonant 
terms v a ^ plus resonant terms v^ b = rjy* a r N *^/(E — Mj^„), where a, b are the considered 
meson-baryon channels, Mjy, is the bare mass of the N* state, and r^*^ describe the iV* — > a 
decays. Such a Hamiltonian can be derived from effective Lagrangians using a unitary trans- 
formation method developed in Ref. j^. By using the two potential formulation j^|, as 
also derived explicitly in Appendix A, one can cast exactly the transition amplitude T a ^(E) 
for the a — > b reaction into a sum of non- resonant t a ^(E) and resonant t^ b (E) terms: 



T a , b (E) = t a , b (E)+t* b (E). (1) 



The first term of Eq. ((TJ is determined only by the non-resonant interactions 

t a ,b{E) = v a , b + Va,c G C {E) t Cjb (E) , (2) 

c 

where G C (E) is the propagator of the meson-baryon state c. The resonant term is 

tU E ) = E rjv*, a (£) [G N *(E)] id r Nhb (E). (3) 

N* .N* 

The resonant amplitude in Eq. (JHJ) is determined by the dressed vertex 

r N *, a (E) = T N *, a + ^N*,b G b (E) t b , a (E) , (4) 

b 

and the dressed iV* propagator 

[G N \Ey%(E) = (E — M° N: )5 i7j - E i7j (E) . (5) 
Here the N* self-energy is defined by 

^(E) = J2r N: ,aG a (E)T^ ha (E). (6) 

a 

In this work, we make the following simplifications. We keep only three channels 7./V, 
KY and irN, and neglect the terms with electromagnetic coupling strenghts higher than 
the first order e. We further assume that the iV* propagator Eq. (JSJ) can be replaced by a 
simple phenomeno logical Breit-Wigner form. Then, Eqs. are reduced to the following 

expressions for calculating the 7iV — > KY and ttN — > KY amplitudes: 

T lN)KY {E) = t jN , K y(E) + E E Jm7* N +^ Y (E)/2 ' (7) 
T nN , KY {E) = Un, ky {E) + E p _ Z'* N , /0 > ( 8 ) 



N* 



E-M N * +iV tot {E)/2' 



with 



rjv* i7 AT — rjy, aN + ^N^yGkY ^N*,KY + tjN,TrNG n NT N *^ N ) , (9) 
Tn*,tvN = Tn*,tvN + [^ttN,KyG KY^ N* ,KY + ^N^nG^nT N * N ] , (10) 

Fn*,ky = r 'n*,ky + [r " n* ,kyGky^ky,ky + 1" 'n*,-knG^js!^n,ky\ ■ (11) 

It is clear that the first step to solve the above equations is to develop models for cal- 
culating all non-resonant amplitudes. To first order in electromagnetic coupling, within the 



considered >yN © KY © tcN space, Eq. (j2J) leads to 



t~,N,KY — V^n^ky[^- + GKY{E)t KY ^KY{E)}+ V 1 N,TvNG 1T N{E)t 1TN ^KY 

= V^n^kyI^- + GKY{E)tKY,KY(E)} + t 1 N,wNG 7r N(E)v 7T N,KY- (12) 

Here we note that the second line of the above equation is obtained from using the well- 
known property vgt = tgv. The non-resonant amplitudes txY,KY and t^^^Y in Eq. (JT2J) 
are obtained by solving Eq. (EI) within the subspace KY © tcN. For numerical reasons, we 
follow the procedure of Ref. to eliminate t W N W N from these coupled equations. We then 
obtain the following equations 



where 



tKY,KY — VKY.KY + ^ V KY,KY Gry tRY,RY, (13) 
KY 

tRY,wN = [VRY,TTN + tRY,RY G RY Vry,ttn] 

X [1 + G^N Kn,ttn], (14) 



V< RY,RY — V KY,RY + ^2 V KY,nN G^N v< ^N,RYi (15) 

ttTV 



with 



,.eff 

■kN.KY 



VttN,KY + X! tnN,nN G w n V^RY- (16) 
ttN 

The pure nN scattering t-matrix t nNt7TN in the above equations is defined by 

tnN,nN = V^ttN + V n N,nN G W N t-nN^N ■ (17) 

We see that Eqs. (fT3j) and (fTTjl are single channel integral equations. The couplings between 
nN and KY channels are isolated in the effective potentials v e ^y KY and v^J KY - Clearly, the 
use of Eqs. (|13| ) - (fT7j) greatly simplifies the numerical task of handling the matrix problem 
associated with the original coupled-channel integral equations in the subspace KY © tcN. 
In fact, this technique will be useful for future investigations including additional channels. 

To solve the above equations, we employ the non-resonant potentials vry,ry, v^n^ry 
derived in Ref. 

.19] 

from effective Lagrangians using a unitary transformation method of 



Ref. [35|. The expressions for these potentials can be found there and will not be repeated 
here. However, we depart from Ref. |19] in two aspects. First, Eq. (17) for determining 







t-wN,-KN was not solved directly in Ref. pj|. Instead, it was estimated from usin g th e empirical 
irN — * 7r iV amplitudes. In this work, we solve Eq. (|T7jl by using v^n^n of Ref. 35| which was 
also derived from effective Lagrangians using the same unitary transformation method. The 
second new aspect of our calculations is to include the distortion effects on the N* decays, 
defined by the term within the square brackets in the right-hand-side of Eqs. ([10 p -([llj h 
which were neglected in the calculations of Ref. . It turns out that these two refinements 
do not change much the quality of the fits to the irN — > KY data. More details will be 
given in the next Section. 

We now discuss the calculation of the non-resonant kaon photoproduction amplitude de- 
fined by Eq. (|12|). While the main contribution to tyN,KY is expected to be from the direct 
transition amplitude v 1 n,ky, the calculations of the coupled- channel effects due to the ttN 
channel require a model for the 7iV — > ttN amplitude t 7 jv )7r jv. The amplitude t lN ^ N is ex- 
pected to be rather complicated in the second and third N* regions. Full construction of 



t-yN,irN is far beyond t 
cal procedure of Ref. 



re scope of this work. To make progress, we follow the phenomenologi- 
4o| to define t 7 /v )7r jv in terms of the empirical — ► nN amplitude and 



27 



28, 



the resonant amplitude constructed from the quark model predictions of Ref. 
Explicitly, we define 

+ _ t-^SAID ,QM,R / 18 \ 

L TrN,-yN — -L ttN~(N l -nN^Ni \ lo ) 

where t^^ N is the quark model amplitude given explicitly in Ref. 42 1 and T^j^P N is obtained 
from the 1995 solution of the SAID analysis. As an alternative, we can replace t^'^ N 
by £^]v 7 jv which is the 7^ — > N* — > tcN amplitude defined by the resonance parameters 
listed by PDG. Unfortunately, the parameters of 7^ — > N* for most of the considered iV* 
are not well determined by PDG. In fact, this work is one of the possible ways to learn 
about these 7iV — > N* amplitudes by considering the photoproduction channels other than 
the ttN channel. We thus use Eq. ()18)) in this work. 

Eq. (fTH|) only defines the on-shell values of the amplitude t^N^N- For the calculation of 
Eq. (fT2"|) . which involves integrations over off-shell matrix elements, we define the following 
off-shell behavior 

t^N^Niq, k , W) = t nNtlN (q , k , W) f} g ' ^ , (19) 

? {Qo, AJ 
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with 

( A 2 \ 2 

F(g ' A) = (a*T?) ' (20) 

_ [(W 2 -m 2 N - ml) 2 - Am 2 ,™ 2 } 1 / 2 
% - ' (21) 

where W is the invariant mass of the tcN system, q is irN off-shell momentum, fc is the 
on-shell momentum of the initial 7 A/" system, and the cutoff A is an adjustable parameter 
in our fit to the -> KY data. We find A = 1.5 GeV. 



B. Direct channel 



For the non-resonant 7 A" — > KY transition amplitude v^n,ky and the resonant amplitude, 



we follow the procedure of Refs. |22|, |2^J . The details can be found there and will not be 
repeated here. Below, we summarize the main points needed in the subsequent Section. 

The chiral constituent quark approach is based on a low energy QCD-inspired La- 
grangian j^j], where the scattering matrix for the photoproduction of pseudoscalar mesons 
can be derived 01 as 

Mfi = (y) |j,.,,|ivj +E {WMgi + mimm^A] 

■ I Ei + lo — Ej Ei — uj m — Ej J 

+M T . (22) 

Here, Ni(Nj) is the initial (final) state of the nucleon, uj(u> m ) represents the energy of incom- 
ing (outgoing) photons, and H m and H e are pseudovector and electromagnetic couplings, 
respectively, and Nj is the intermediate baryon. 

The first term in Eq. IS du SCB. gull term. The second and third terms correspond to 
the s- and ^-channels, respectively. The last term .Mr is the t-channel contribution. 

The contribution from the s-channel resonances to the transition matrix elements can be 
written as 

9 A/f k' 2 +q 2 

kaCQM _ ZMn* -^r- . , , 

Mn * ~ W*-M N .(M N .-iT(q)) 6 AN *> (23) 

with k = \k\ and q = \q\ the momenta of the incoming photon and the outgoing meson, 
respectively; W is the total energy of the system, e~( h2+q2 ^ 6a h° a form factor in the harmonic 
oscillator basis with the parameter o!^ lo related to the harmonic oscillator strength in the 

8 



wave-function, and Mn* and T(q) the mass and the total width of the resonance, respectively. 
The amplitudes An* are divided into two parts: the contribution from each resonance below 
2 GeV (these transition amplitudes have been translated into the standard CGLN amplitudes 
in the harmonic oscillator basis), and the contributions from the resonances above 2 GeV , 

n 

which are treated as degenerate f51j. 



22 1 a new set of real 



The contributions from each resonance is determined by introducing 
parameters Cn* for the amplitudes An*'- 

An*^C n *A n *i (24) 

so that 

Mn* = C 2 N ,M C N 9 M , (25) 

where M. e £* is the experimental value of the observable, and A4 N 9 M is calculated in the quark 
model |23j. For instance, for resonance with mass < 2 GeV, the SU(6) ® 0(3) symmetry 
predicts Cn* = 0.0 for Sii(1650), -Di3(1700), and Dis(1675) resonances, and Cn* = 1-0 for 
other ones. However, deviations from those central values are anticipated within the broken 
SU(6) <S> 0(3) symmetry, due for example to one-gluon exchange mechanisms 



III. RESULTS AND DISCUSSION 



This Section is devoted to the application of our formalism to the production of kaon- 
hyperon final states in txN and 7p collisions. 

To that end, we need first to study ttN — > KY, and KY —>■ KY processes. In the 
following we first compare our ttN — > KY results with the relevant data and also extract iV* 
information within the considered model. Then we present results for the photoproduction 
channel and discuss issues related to the missing resonances. 



A. irN -> KY Reaction 



As seen in Eq. (|12j) . to calculate 7iV — > KY amplitude our first step is to construct 
the non-resonant amplitudes txY,KY and t n N,KY- These are obtained within our model 
by solving the coupled-channel equations (fT3 ^) - (fT7j) . The input of these coupled-channel 
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Notation 


Resonance 


Coupling 


Value 






Ik an 


-0.61 






/kejv 


0.12 






/ttEA 


0.08 






/vrEE 


0.00 


A'i 


5u(1650) 1/2 - 


/if AJV4 


-0.25 






/lfEiV4 


-0.20 


iV5 


£>i 3 (1700) 3/2" 


/KAJV5 


-0.33 






}kT,N5 


0.08 


N6 


Pn(1710) 1/2+ 


IkAN6 


0.09 






}kT,N6 


-0.32 


N7 


Pi 3 (1720) 3/2+ 


Skani 


-0.56 






}KT,N7 


0.54 


Dl 


#31 (1900) 1/2- 


/jfSDl 


0.09 


D2 


P 3 i(1910) 1/2+ 


/jfSD2 


0.20 


D3 


P 33 (1920) 3/2+ 


}kt.dz 


-0.20 


L3 


Sox (1670) 1/2- 


fnSL3 


-0.20 


Lb 


Poi(1810) 1/2+ 


/vrEL5 


-0.01 


SI 


Pn(1660) 1/2+ 


/ttASI 


-0.20 






/ttESI 


-0.20 


54 


Di 3 (1670) 3/2- 


/ttAS4 


0.22 






/ttES4 


0.05 



K*NY couplings 


fV 

Jk*na 


0.71 




Jk*na 


-3.97 




fV 


-0.53 




J if* ATE 


0.52 


cut-offs 


As 


623.0 




A u 


1468.0 




A< 


930.0 




A W N 


1491.0 


off-shell 


X 


2.0 


Reduced \ 2 




1.86 



TABLE I: Coupling constants in ttN — > KY and KY — > KY . The values are extracted from our 
minimization procedure. The parameters are defined in the model B of Ref. 



equations are the potentials vky,-kN, vky,ky, and an effective non-resonant amplitude t^N^N 
which is defined by Eq. ()17j) . The parameters of these potentials are then adjusted along 



with the N* parameters associated wit 



and rr p —>■ K°H° data 



53, 



54 



55 



i the resonant term of Eq. (|SJ) to fit the rr p — > K°A 



This policy was pursued in Ref. but with the simplifications that the distortion fac- 
tors, the terms within the square brackets in Eqs. (|TU ]) -([TT |) . were not included in calculating 
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the resonant term of Eq. (jHJ). Furthermore, the non- resonant i^N^N defined by Eq. (fTTjl was 
only roughly estimated using the empirical ttN amplitude. 



I ' 1 ' 

1 W= 1.661 GeV 

— t It 


It 5' „ 

W= 1.683 GeV ft 1 
F 1 1 1 1 1 1 1 


1 | 1 | 1 | £1 

W=1.724GeV <£ 
5 1 1 1 1 1 1 1 


' 1 ' 1 ' 1 

_ W= 1.758 GeV ^ 



_ W= 1.847 GeV _ 

,1,1 


- W=1.879GeV -r 


1 W=2.059 GeV 

• Saxon 

O Baker 

* Knasel 


- W=2.159GeV 


_ W=2.259 GeV I 



'-1 -0.5 0.5 -1 -0.5 0.5 -1 -0.5 0.5 1 



cos(6) 



FIG. 1: Differential cross-section for the reaction 7r p-> K°A. The solid curves are from the fits 
using the coupled-channel model of this work. Data are from Refs. J53, 55]. 



In this work, we have corrected these two deficiencies as discussed in Section II. A and 
thus have refined the potentials tV7v,ify and vky,ky and the relevant N* parameters. 

The fitting procedure is explained in detail in Section III of Ref. |l9j]. Here we recall 
a few points to make the present Section self-consistent. In that paper, we classified the 
parameters in three sets (Tables I to III in Ref. Q]). Set I includes 9 couplings, the values of 
which are taken from the SU(3)-symmetry predictions or PDG partial decay widths; namely, 
Unn, Unn*,Una*, and / ffjvs ., with N* = S n (i650), D 13 (1700), Pn(1710), Pi 3 (1720), 
A* = 501(1670), P i(1810), and £* = P n (1660), D a3 (1670). 

The adjustable parameters are in the remaining two sets. Set II includes the following 
coupling constants: /kyn, Ikyn*, Iky a*, fnYY, and f n YY*- The values extracted for those 
parameters within the present work are given in Table HJ rows 2 to 22. Here we follow model 
B of Ref. [3] by allowin g th e parameters of the model to v ary by ±20% around the central 
values taken from PDG |36| and/or from quark model 27|,|28j predictions. Finally rows 23 
to 31 in Table |l] correspond to Set III in Ref. Q|. 

In Figs. ^to|3J the results of our model are compared with the differential cross-section 
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cos(6) 



FIG. 2: Differential cross-section for the reaction tt p ^ The solid curves are from the fits 

using the coupled-channel model of this work. Data are from Refs. I54I 



53 



54. 



56 



571 ] for processes it p — > K°A and 



and recoil hyperon polarization data 
n- p -> K°T.°. 

In Figs. ^ and 121 we show the quality of our fits to the differential cross section data for 
ix~p — > K°A and n^p — > fT°£ , respectively. In Figs. El and 0] our results for the asymmetry 
data for the same reactions are depicted. The acceptable agreement between model and 
data, as we., as , L „ compare we., with our previous resu.ts Q. Nevertheless, we consider 
the present model slightly more reliable than the model B in Ref. [19]. Actually, some of 
the coupling constants, Table HI get (much) closer to constituent quark model values [29] , 

e-g- /jCADi 3 (1700)) /i<"SI?i3(1700)) /fir APu (1710) > /fi. AP 13 (1720) > /iCS5 3 i (1900) , and /fi-SP 33 (1920)- 

In summary, it turns out that the aforementioned improvements do not change much with 
respect to our previous model Q]. For the KY — > KY processes, we also get comparable 
results to those reported in the latter paper. There is no data for KY — > KY scattering 
to test our model. The situation might change in the near future with the advent of highly 
accurate data from the EPECUR 58 1 and J-PARC 3 

projects. Those data will certainly 
afford deeper insights into the meson-baryon interactions. However, the results shown in 
Figs. HUH are sufficient for the purpose of studying coupled-channel effects. 
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< 




0.5 0.5 1 



FIG. 3: A recoil polarization asymmetries for the reaction tt p — ► K°A. The so 
the fits u: 



id curves are from 




0.5 o 0.5 1 



FIG. 4: £ recoil polarization asymmetries for the reaction ir p — > K°T,°. The solid curves are from 
the fits using the coupled-channel model of this work. Data are from Ref. |57| |. 



B. 7p — > K + A Reaction 

We have performed a thorough study of all the latest relevant data fTablelTT|). The data 
released in December 2005 by LEPS p are not included in our fitted data base. However, 
they are depicted in the relevant figures below and briefly discussed. 
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The strong interaction channels amplitudes v w n,ky and txy,KY are determined above, 
and t^N,nN computed from Eq. (|I5jl. 

Lab/Collaboration Observable # of data points Ref. 

ELSA/SAPHIR Differential cross-section 720 [3j 

JLab/CLAS Differential cross-section 1068 [1] 

JLab/CLAS Recoil polarization asymmetry 233 [2] 

SPring-8/LEPS Polarized beam asymmetry 44 [5] 

Bonn synchrotron Polarized target asymmetry 3 [6bJ 

TABLE II: Data sets investigated in the present work. Here, we have not included 268 cross-section 
data points from CLAS for E 1 > 2.6 GeV (W > 2.4 GeV), in order to concentrate on the baryon 
resonances energy range. 



For 
model 



3oth non-resonant and resonant 7^ — > K + A amplitudes we use a constituent quark 
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601 ]. We recall that the resonant term of Eq. (J7J) contains a term 



tyN,KY - V jN,Ky[^ + G KY t K Y,KY], (26) 

with 

R _ r|v*, 7 iv T N . iKY 

V lN ,KY~ E _ MN . +ir tot( E)/2 - 

To use the N* contributions defined by Eqs. (J2HJ)- (J23J), we replace the above expression by 

V ^N KY = Cjy*Mx? M , (28) 



where M^? M is calculated 



231 ] from the constituent quark model. The SU(6) (g> 0(3) 



symmetry breaking coefficient Cn*, Eq. (|2B|) . are treated as constrained adjustable param- 



eters 



22, 



24j in fitting the data. 



1. Model search 

In this Section, we explain the procedure used to build a model for all available data. 
Here we would like to emphasize that the CLAS [2] and SAPHIR [IJ data released in 2004, 
with some 2000 data points for differential cross-sections, showed significant discrepancies 
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with each other. This fact led the phenomenologists either to concentrate on one of the two 
sets or produce one model per data set. This uncomfortable situation is now significantly 
cured thanks to the CLAS Collaboration's new data |l|, made available in 2005. In an earlier 
attempt j2l|, we underlined this improvement in experimental data base and reported our 
preliminary results obtained with respect to both SAPHIR 2004 and CLAS 2005 data. 

Table ITT1 summarizes the content of the data base used to determine the adjustable param- 
eters of our approach; namely, known resonances strengths. Additional parameters due to 
the introduction of new resonances are discussed in Section liHB 21 Differential cross-section 
data provide, of course, the main constraints on the model ingredients. Consequently, our 
starting point was to fit separately the CLAS and SAPHIR cross-section data, for which the 
reduced x 2 s turned out to be 2.1 and 1.3, respectively. The significantly larger Xd.o.f f° un d 
using the CLAS data is due to their smaller uncertainties compared to those of SAPHIR 
data. However, this fact might not be the only source of the difference in % 2 s. Actually, two 
considerations are in order here: 

i) The earlier data from CLAS |2] showed significant discrepancies with SAPHIR jsj data. 
Although the new CLAS data set has significantly reduced those discrepancies, in some 
phase space regions results from the two data set differ still by more than 2cx; 

ii) The strengths of resonances, which constitute our main adjustable parameters, are 
rather tightly constrained by SU(G)®0(3) symmetry. Consequently, the fact that we obtain 
a much better Xdo f f° r the SAPHIR data compared to the one for the CLAS data leads to 
the conclusion that our approach is more in line with the SAPHIR differential cross- sect ion 
data than with CLAS results. 

Keeping the above considerations in mind, we present two models here: 

i) Model M\. all SAPHIR and most recent CLAS differential cross-sections (first two rows 
in Table II) were fitted simultaneously. 

ii) Model M2. all cross-section and polarization asymmetries (Table ITT]) were fitted simulta- 
neously. 

Extracted values for the eleven adjustable parameters are given in Table ITTT1 That Table 
contains the KYN coupling constant and the strengths of known resonances with masses < 2 
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GeV. The higher-mass, known resonances are treated as degenerate in a compact way 
and bear no symmetry breaking coefficients. Moreover, the Roper resonance, although 
explicitly present in our approach, does not contribute to the reaction mechanism due to its 
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low mass with respect to the reaction threshold. In addition to those known resonances, we 
also introduce 3 adjustable parameters per each of three newly proposed Su, P13, and Di 3 
resonances, as discussed below (see Table HV1 in Section llll B 2|) . 



Parameter Model M x Model M 2 



9KNA 


8.02 


8.00 


Csii(1535) 


-0.85 


-0.82 


C5n(1650) 


-0.10 


-0.22 


Cpn(mo) 


1.79 


-1.08 


Cdi 3 (1520) 


-2.00 


-2.00 


Cdi 3 (1700) 


0.16 


-0.19 


Cpi 3 (1720) 


-0.40 


0.05 


Cpi 3 (1900) 


0.80 


1.60 


Cdi 5 (1675) 


-0.09 


0.22 


Cfi 5 (1680) 


1.43 


1.99 


CFi 5 (2000) 


1.28 


1.59 


Xd.o.f 


2.49 


3.32 



TABLE III: Kaon-nucleon-hyperon coupling constant, SU (6) (g>0(3) symmetry breaking coefficient 
Cn* as in Eq. ()28|) . and reduced % 2 for models Mi and M 2 . 

Here we would like to comment about the extracted values of adjustable parameters. 
The coupling constant gxNA is very close to its lowest limit within broken SU(3)- 
symmetry. This parameter, like several other adjustable ones, is driven by CLAS data. 
Actually fitting only the SAPHIR data leads to qk n a =9. 70. Finally, the Xd.o.f f° r the model 
Mi is significantly higher than obtained by fitting only the SAPHIR data. Actually the 
integrated % 2 for the latter data set increases by more than a factor of 2, i.e. the adjustable 
parameters are driven by the CLAS data. However, in going from the model Mi to M2 
that integrated x 2 stays stable, while the integrated x 2 f° r CLAS data increases by roughly 
30%. Moreover, in the integrated x 2 s for the models Mi and M 2 , CLAS data represents 
roughly 55% and 48%, respectively, while SAPHIR data account for about 45% and 29%, 
respectively. These results indicate that, within our approach, the SAPHIR data show larger 
compatibilities with the polarization data, than does the CLAS data. 
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FIG. 5: Differential cross-section for the reaction 7p — > K + A as a function of total center-of- 
mass energy. Dotted and solid curves correspond to models Mi and Mo respectively. Data are 

li Q Q 

from Ref. [1[ (open diamonds), Ref. [3|] (full circles), and LEPS la] (open squares in the cells 
corresponding to cos# = 0.75 and 0.85). Plotted data from Ref. [l| are measured at cos(9)' = 
cos{9) + 0.05. 

In Figs El to |H1 results for models M\ and M2 are compared with the most recent data. 

In Fig. El excitation functions at 19 angles, for 9% ~ 25° to 150° are shown as a function 
of total center-of-mass energy for W = 1.6 GeV to 2.3 GeV. Except in very few phase space 
regions, the two models give identical results. Given the discrepancies between the two fitted 
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cos(9) 



FIG. 6: Angular distribution for recoil A polarization asymmetry for the reaction *yp — ► K + A. 
Curves as in Fig. 0D Data are from Ref . [3] • 

data sets, our models give an acceptable account of the differential cross-sections. In the 
same figure, we show also the very recent LEPS data [f| for cos# = 0.75 and 0.85. They 
turn out to be closer to the CLAS data, rather than to SAPHIR results. 

With respect to the polarization observables, we recall that model Mi ( dotted curve) has 
been obtained by fitting only the cross section-data. So, in Figs El to El dotted curves are 
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predictions. While the full curves (model M2) result from fits to differential cross-section 
and polarization observables data. 
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FIG. 7: Angular distribution for polarized beam asymmetry for the reaction 7p — > K + A. Curves 



as in Fig. OD Data are from Ref. 



(circles) and Ref. p (squares). 



In Fig. |H1 angular distribution of polarized recoil A asymmetry is depicted for W ~ 1.6 
GeV to 2.3 GeV. Models Mi and M 2 give significantly different results and the latter model 
reproduces the data quite well, except for a few lowest energy ones. It is worthwhile mention- 
ing that although recoil data represents less than 10% of the data base points, and contribute 
to the total x 2 by the same percentage, they have a significant effect in the determination 
of the model ingredients. 

The polarized photon beam asymmetry, Fig. [7[ data stand for less than 2% of data 
base points, but generate about 13% of the total \ 2 - We recall that the fitted data come 
from Ref. p| and are shown as open circles in all 9 cells of Fig. while the very recent 
data depicted as open squares, were not included in the fitted data base. Here, model 
M 2 (solid curves) shows an improvement over Mi (dotted curves) when compared with the 
data. According to our results, further measurements of this observable around 9k ~ 90° 
would put strong constraints on the models search. 



Polarized target asymmetry has been measured only by one group 



about 3 decades 
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ago. For completeness, we compare our models with those few data points, Fig. [SI showing 
that the model Mi gives a better agreement with those data. Contribution of those data to 
the total x 2 is around 0.1%. 




FIG. 8: Excitation function for Polarized target asymmetry for the reaction jp —* K + K. Curves 
as in Fig. OD Data are from Ref. [Q2 1 . 



In summary, the model M2 provides a reasonable description of the whole data base 
Figs. [SHBl and the comparisons of the resulting parameters listed in Table III indicate the 
importance of having polarization observables data in the study of N* resonances. 

2. Search for new resonances 



For about three decades, several approaches have been predicting baryon resonances not 
seen in extensively investigated 7rN channels, issues related to those missing resonances 
have recently been reviewed 2(3, [ill Q 

The search for missing resonances has been initiated by predictions formulated in three pi- 
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3,ii 



ii) Algebraic approach [3l| , 



oneer approaches: i) Relativized quark formalism 
iii) Hypercentral constituent quark model (3^. 

Moreover, several authors have reported about three missing resonances with masses 
between 1.8 and 2 GeV, namely, Six, P13, and D 13 . In the present work, we have investigated 
possible contributions from such resonances to the jp — > K + A reaction mechanism. Before 
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presenting our results, we give a brief account of findings reported in the literature with 
respect to those resonances. It is worthwhile keeping in mind that all the results mentioned 
below and referring to the CLAS data, use the CLAS 2004 results 2] and not the more 
recent ones jl|. So, conclusions based on those works have to be updated, 
a) Third S n 

The extracted values for the mass and width of a new Sn are close to those predicted by 
the authors of Ref. [39] (M = 1.712 GeV and T = 184 MeV), as a KY bound state. 



40] 



23 M in the 



The chiral constituent quark approach used in the present work served 
interpretation of the 7p — > r\p data and put forward strong indications for a third Sn with 
M = 1.780 GeV and T=280 MeV. 



For the one star 5n(2090) resonance [36j and where the mass ranges between 1.880 GeV 



4l| produces the following values 



to 2.180 GeV, the Zagreb group's coupled channel analysis 
M = 1.792 ± 0.023 GeV and T = 360 ± 49 MeV. The same one star resonance was invoked 
in the 1.932 GeV to 1.959 GeV range, using a reggeized isobar model to investigate the 
7P — > rfp reaction. Still another isobar approach Q] investigation of the 7p — > r\p puts 
forward an Sn resonance with M =1.825 GeV and T=160 MeV. 

A self-consistent analysis of pion scattering and photoproduction within a coupled channel 
formalism, indicates [3] the existence of a third S\\ resonance with M =1.803 ± 0.007 GeV. 

Finally, one of the main recent experimental sources on baryon resonances comes from 
the BES Collaboration j^l^l, using J/\l/ decay channels. In an early stage, they concen- 
trated (4^ on neutral pion and rj final states: J/^/ — > ppn°, ppr/. The authors could identify 
the two known Sn resonances and extracted their masses and widths in agreement with the 
PDG values. They found a structure at M = 1800 MeV, the quantum numbers of which 
could not be identified because of lack of statistics. 

b) Third P 13 

Very recently, the BES Collaboration has released [3] data for charged pion final states: 
J/ty — > pn + n, nir~p. Besides again identifying the two known Sn resonances, they put 
forward the following interesting results: i) The Roper Pn resonance's mass and width are 
reported, M = 1358 ± 6 ± 16 MeV and T = 179 ± 26 ± 50 MeV, to be significantly smaller 
than their widely used values, ii) A fourth resonance was identified by the authors with M 
= 2068 ± 3l^ MeV and T = 165 ± 14 ± 40 MeV, 3/2+ spin parity. 

c) Third D n 
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The first indication of a new D\z with a mass close to 



and Bennhold 



who interpreted the SAPHIR 1998 data 



.9GeV was suggested by Mart 
within an isobar approach. 



Subsequently, it was shown that those data could be reproduced both within an isobar 
nrodel 0, enrbodyrng oft-sheU effects, and a constituent quark approach Q. Moreover, 
recent data [3( released in 2004 by the SAPHIR Collaboration did not confirm the structure 



reported in their 1998 paper. Afterwards, Mart et al. |66( reached the conclusion that the 
manifestations of such a resonance appeared to be poorly-determined. 

Within an isobar model, including s- and i-channel contributions in the tree approxima- 
tion, Anisovich et al. analyzed the processes — > ttN, r/N,K + A, K + T,°, K°T, + and 
suggested a new D 13 with M = 1875 ± 25 and T = 80 ± 20. The authors report a less strong 
indication for an additional -D13 with M = 2166 and T = 300 ± 65, that they attribute 
to the N*(2080) of PDG. However, recent results from the CB-ELSA Collaboration 3 on 
the 7£> — > N*(A*) — > ir°p puts this latter two star resonance at M = 1943 ± 17 and T = 
82 ± 20. 

A hybrid isobar plus Regge model has been developed by Corthals et al. j3| . According 
to the Regge background model used, a _Di3(1895) appears or vanishes. The authors suspect 
a role for significant final state interactions not included in their approach. Such effects are 
also absent in all isobar models discussed above. 

Such effects, as well as intermediate state reactions, are of course embodied in the coupled- 
channel approaches based on the K-matrix formalism developed by the Giesssen |16| and 
Groningen jl^l groups, though both groups use isobar models for the direct processes. Nei- 
ther of those works show evidences for new resonances. However, the Giessen group fitted 
separately SAPHIR and CLAS 2004 data and the Groningen group used only SAPHIR data. 

Finally, an investigation j47| of the relations between the S'-matrix and time delay in ttN 
interactions, concluded that a -Di 3 (1940) could appear. 

Given the results from other investigations outlined above, we proceed to the presentation 
of our findings with respect to those new resonances. 

In Section riHB 1| we presented our model and made comparisons with all available data 
sets. In this Section, we use the model M 2 discussed in Section lHI B 11 in order to investigate 
possible manifestations of three missing resonances: Su, P13, and Di 3 . For that purpose, 
we have attributed 3 adjustable parameters (mass, width, and strength) to each of those 
resonances in the minimization procedure. The extracted parameters are given in Table ITVl 
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New resonance 


Property 


Model Mi 


Model M 2 




Mass 


1.833 


1.806 


Sn 


Width 


0.288 


0.300 




Strength 


0.40 


0.15 




Mass 


1.974 


1.893 


Pl3 


Width 


0.108 


0.204 




Strength 


0.12 


0.28 




Mass 


1.912 


1.954 




Width 


0.316 


0.249 




Strength 


1.50 


0.98 



TABLE IV: Determined parameters for the third Sn, P13, and D\s resonances. 

In order to ascertain the role played by each additional resonance given in Table IIVI 
we proceed as follows. In Figs. E] to ^] we show the same observables as in Figs. |3] to El 
respectively. For each observable, the model M 2 is depicted again. The three other curves 
in Figs. El t oE3 correspond to the model M 2 with one of the additional resonances switched 
off, without further minimizations. In those figures, the curves are: M 2 without the third 
Sn (dotted curve), P13 (dot-dashed curve), and Di 3 (dashed curve). 

From the differential cross-sections (Fig. EJ) we infer that the 3 rd Sn has a significant role 
in the backward hemisphere and the effect gets enhanced in going to most backward angles. 
The manifestations of this resonance vanish for W < 1.9 GeV. Moreover, the interference 
terms due to this resonance appear to be destructive in the full model M 2 . 

Contributions from the 3 rd P13 resonance are confined roughly to the energy range 1.8 < 
W < 2.0 GeV with increasing magnitude in going from forward to backward angles. Those 
contributions are rather small, but non vanishing in the whole phase space. 

The most significant effects due to the 3 rd D13 resonance are around $k ~ 90° and 
W ~ 1.9 GeV. The interference terms come out to be constructive in the forward hemisphere 
in the whole energy range and in the backward hemisphere for roughly W < 2.0 GeV. 

The recoil hyperon polarization asymmetry, Fig. E3 shows no significant sensitivity to 
the third Sn and P13 except in very limited phase space regions, while switching off the 3 rd 
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W (GeV) W (GeV) 

FIG. 9: Differential cross-section for the reaction ~fp — > K + A as a function of total-center-of-mass 
energy. Solid curve corresponds to the full model Mi. Dotted, dot-dashed, and dashed curves 
correspond to the full model without the 3 rrf Sii, 3 rrf Pi3, and 3 rrf -Di3, respectively. Data are as in 
Fig. El 

Di3 leads to important variations in the model values for roughly W > 1.9 GeV, mainly in 
the forward hemisphere. 

The same trends are observed for the polarized beam asymmetry with respect to the third 
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FIG. 10: Angular distribution for recoil A polarization asymmetry for the reaction — > K + A. 
Curves as in fig. El Data are from Ref. Q|. 

Six, Fig. ^2 The highest sensitivities to the two other resonances appear in the backward 
hemisphere and are significant for the 3 rd -D13. 

For the sake of completeness, in Fig. ^] we show the excitation function at Qk — 90° for 
the polarized beam asymmetry. As already mentioned, this observable is by far the least 
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FIG. 11: Angular distribution for polarized beam asymmetry for the reaction ~Jp — > K + A. Curves 
as in fig. ED Data are from Ref. 

studied experimentally. Our results might nevertheless indicate that the 3 rd -D13 produces a 
significant structure at higher energies. 
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FIG. 12: Excitation function for polarized target asymmetry for the reaction jp — > K + A. Curves 
as in fig-IU Data are from Ref. [61 1. 
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3. Role of resonances in total and differential cross-section, and polarization observables 



Total cross-sections have been extracted by both CLAS and SAPHIR collaborations. 
Those data were not included in our fitted data base. The postdiction of our model M 2 is 
depicted in Fig. EH in bold full curves. In each of the four cells, we show in addition the 
results of that model M 2 with only one resonance switched off at a time. 
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FIG. 13: Total cross-section for jp — > K + A as a function of total center-of-mass energy. The 
bold full curves come from the model M%. In each cell different curves correspond to the model 
M2 with one resonance switched off, as singled out in each cell. Data are from Refs. Q, j|. 

The first observation concerns the discrepancies between the two data sets. As already 



pointed out [2l|, here the discrepancies are more significant than in the case of differential 
cross-sections. This increased discrepancy is likely due to two facts: i) the two collabora- 
tions have performed measurements in non completely overlapping phase space regions, ii) 
different extrapolation methods to the unmeasured angular areas are used. The total cross- 
section extracted from differential cross-sections might then be misleading if included in the 
data base and/or used to draw strong conclusions about the reaction mechanism. There is 
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however a puzzling point, namely, total cross-sections for the 7p — > K + Y,° extracted by the 
same collaborations, agree quite well with each other (e.g. see Figs. 20 and 21 in Ref. 

The model M 2 ingredients are dominated by both data sets up to W ~ 1.7 GeV, by 
CLAS data up to W ~ 2.0 GeV, and by SAPHIR data above that region. Two structures 
appear at about 1.7 GeV and 1.9 GeV. 

To gain better insight into the role played by each resonance with mass M < 2 GeV, we 
show curves obtained using the model M 2 by switching off each resonance. Table IS gives 
the x 2 f° r each case, without further minimizations. 

In the following, we concentrate on the model M 2 in order to investigate contributions 
from various resonances. The points discussed below do not depend on the total cross-section 
data, but they embody effects from all other fitted observables. Moreover, we present the 
effects of each resonance with respect to the fitted data base. Here, in order to limit the 
number of figures, we summarize our findings in Table 13 the content of which is explained 
below. 

The integrated \ 2 m model M 2 can be written as a sum of five partial xi s i 

X 2 = EX, 2 , (29) 

i=l 

where i refers to the data sets, namely, 

= 1 : CLAS differential cross-sections, (da)cLAs', 
= 2 : SAPHIR differential cross-sections, (d<r ) saphir) 
= 3 : CLAS recoil polarization asymmetry, P; 
= 4 : LEPS polarized beam asymmetry, S; 
= 5 : Bonn polarized target asymmetry, T. 
Then, for each switched off resonance, and without further minimizations, we obtain the 
relevant integrated [x 2 ]m 2 -n* and partial [xJ\m 2 -n* for the observable i. (Here, the subscript 
M 2 — N* denotes that the particular resonance N* has been turned off.) Finally, we define 
the following ratio: 



[Xi\M 2 

which gives a measure of the role of the relevant N* with respect to the observable numbered 
i. In columns 3 to 7 of Table IS following found intervals are reported: 
1.0 < Ki < 1.5 : * 
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1.5 < Ki < 1.8 : ** 
2.0 < Hi < 4.4 : *** 
6.0 < Ki < 8.0 : **** 
10.0 < Hi < 13.4 : ***** 

Thus more stars indicate a larger role for a particular resonance on a particular observable 

i. 



Switched off N* (xM 2 -N*)d.o.f. K {da)cLAS n {da)sAPHIR TZ P Tlx TZ T 



<bii(15oo) 


1 n o 

1U.3 


^ >(< 


^ >(< ^ 






* 


Dii(lb5Uj 


5.7 


H< ^ ^ 




* 




* 


511(1806) 


6.5 


*** 


*** 


** 


* 




Pll(1710) 


3.3 


* 


* 


* 




* 


Pis(1720) 


3.4 




* 


* 


* 




As (1900) 


13.0 


**** 




*** 


*** 




Pis(1893) 


4.6 


** 




* 






Di 3 (1520) 


20.0 




*** 


*** 




** 


^13(1700) 


3.5 


* 


* 


* 




* 


£>is(1954) 


26.4 




*** 


*** 


*** 


* 


£>is(1675) 


3.6 


* 


* 




* 


* 


Fi 5 (1680) 


5.0 


** 


* 


** 


* 




F 15 (2000) 


7.1 


*** 


* 


* 


*** 


* 



TABLE V: Schematic presentation of the role played by each resonance in the process 7p — ► K + A. 
First column: switched off resonance in model M2; second column: reduced x 2 without further 
minimizations to be compared with the (xM 2 )d.o.f = 3.3 for the model M2 (see Table ITTT|) . The 
third to sevenths columns give the intervals of TZi (Eq. I30|) with the number of stars as defined in 
the text. The three new resonances investigated here are given in bold. 

In few cases, the TZi is slightly smaller than 1.01, shown by a hyphen (-) in that Table. 

The cell on left-top (Fig. IT3*|) shows the effects of Su resonances. The lightest resonance 
affects the total cross section significantly above its mass, due to constructive interference 
terms, and contributes clearly to the first maximum. This is also the case for the Su (1650), 
with smaller effects close to threshold. The 3 rd Su intervenes around 1.8 GeV and brings in 
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destructive interference. The first and third Sn resonances play important roles f Table IV)) in 
the differential cross-section data from CLAS and SAPHIR, while the second one is present 
only in the CLAS data. 

For the P-waves (Fig. 1131 left-bottom cell), Pn(1710) and Pi3(1720) have negligible con- 
tributions and they do not appear in any of the observables (Table IVj). According to the 
same Figure and Table, the Pi3(1900) has strong manifestations within the CLAS differential 
cross-sections and, to a less extent, in the P and T polarization observables. 

In the spin 3/2 D- waves case (Fig. right-top cell), the first such state plays an impor- 
tant role with interference effects turning from constructive to destructive around 1.9 GeV. 
Table IVI shows that the Pi 3 (1520) is a crucial ingredient in reproducing the CLAS data and 
is important with respect to the SAPHIR results. The role of the Pi3(1700) is negligible, 
while the 3 rd P13 has a clear role between roughly 1.8 GeV and 2.0 GeV (Fig. IT3"j) and turns 
out to be a key element, Table fVl for all observables, except T. 

The spin 5/2 D- and F- waves show no significant effects in the total cross-section (Fig. El 
right-bottom cell). However, Table underlines the importance of the P15 resonances, 
especially the second one. 

To summarize this Subsection, we find that: 

• Among the known resonances, the most relevant ones are: 5*11(1535), Pi3(1900), and 
#13(1520); 

• Three other ones are required by data other than those from SAPHIR: S , n(1650), 
P 15 (1680), and P 15 (2000); 

• Among the three new resonances, the Pi3(1954) plays a crucial role in all observables, 
except perhaps in the beam polarization asymmetry. The 511(1806) plays an important 
role with respect to both differential cross-section data sets, and the polarized recoil 
data. The Pi3(1893) has a less strong role than the two previous resonances and shows 
up mainly in the CLAS cross-section data. 

As mentioned above, all the curves with one resonance removed and depicted in Fig. 
are obtained without further minimizations. For the sake of completeness, we shortly report 
about all possible configurations, after minimization, without new resonances, with one or 
two of them included. 
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The reduced % 2 s are given in Table IVII The second column shows the result with only 
known resonances, with xio./ = 5-7. Adding either a third Sn or P i3 decreases the Xd.o.f 
by roughly 18% (3 rd and 4 th columns). While a third -D 13 improves the x\of D Y about 
32%. Columns 6 to 8 show the effects of combinations of two new resonances. The smallest 
Xdof i s obtained by the Su and Di 3 pairs. The last columns recalls the result for model 
M.2- It is interesting to mention that the mass of those resonances stay stable within 50 
MeV through the 7 configurations. These results support our conclusions above on the a) 
important role played by the third D 13 , b) improvement due to an additional Sn, c) less 
significant contribution from P13. 



New resonances — 


* None 


S 


V 


V 


SP 


VD 


SD 


SPD 


Xd.o.f 


5.7 


4.8 


4.7 


3.9 


4.2 


3.6 


3.4 


3.3 



TABLE VI: Dependence of the x\ Q f on a ^ possible combinations with respect to the three new 
resonances. 

Finally, we outline here the results obtained by using only the direct-channel calculation 
(with all multi-step processes turned-off). Embodying only the known resonances leads to 
Xd.o.f ~12, to be compared to 5.7 in Table IVII The Xd.o.f S e ^ s improved by adding new 
resonances and goes down to ~4 when the 3 new resonances are included. Although the 
Xd.o.f § e ^ s an acceptable value, some of the adjustable parameters turn out to be irrealistic. 



4- Coupled- channel effects 

It is important to illustrate here the differences between the coupled-channel ap- 
proach presented here and the often used approximations in the literature: i) the tree- 
diagrammodels (direct-channel) neglecting multi-step phenomena, ii) the coupled- channel 
K-matrix approaches neglecting off-shell effects. 

A tree-diagram model can be obtained from the formulation presented in section [H] by 
turning off all multi-step processes. Namely, the tree-diagram amplitude is simply 

TyN^KY = VyN,KY + V^ N>K y , (31) 

where v^n,ky is the non-resonant amplitude and v^ N KY ls the resonant amplitude calculated 
from Eq. (JUJ). 
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FIG. 14: Differential cross-section excitation functions at four angles for jp — > K + A. The curves 
are: model M2 (full curves), direct-channel results obtained by turning-off multi-step processes in 
the full calculation (dotted curves), and off-shell effects swithed-off in the full calculation (dashed 
curves). Data as in Fig. [5J 

The importance of the coupled- channel effects can be seen by comparing the results from 
Eq. (J3*T]) and the coupled-channel equations Eqs. (ffj)- (fT7j) . We see in Fig. H31 that when the 
coupled-channel effects are turned off, the resulting differential cross sections (dotted curve) 
would largely overestimate the cross sections; especially in the energy region W ~ 1.6 — 2 
GeV. Obviously, the resonance parameters extracted from using the tree-diagram model will 
contain such theoretical uncertainties. 

Moreover, within the coupled-channel formalism, the role played by off-shell effects is 
depicted in Fig. El The dashed curves there show our results when the off-shell treatment 
is turned off. Sizeable effects are present in the same energy range as above. 



5. Meson cloud effects on N* excitations 



In the dynamical study of the A resonance, it was found |35|,|38( that the dressed — > A 
transition contains a large contribution due to the mechanism that the bare A state is not 
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directly excited by the incident photon, but by the pion first produced by the non-resonant 
mechanism. This contribution, commonly termed as the "meson cloud effect", can also 
be identified within the coupled channel model considered here. Within the formulation 
presented in Section II, the meson cloud effect is contained in the terms within the square 
brackets of Eq. JSJ). Obviously such a meson cloud effect is absent in the tree-diagram model 
defined by Eq. (p?T]) . We also note that the calculation of these meson cloud terms involve 
integrations over the off-shell matrix elements of non-resonant amplitudes t 7 7v,/£v and t^N^N- 
Such off-shell dynamics is neglected in the K- matrix coupled-channel model jlai ]. 




FIG. 15: The multipole amplitudes calculated from the — > N* — > KY resonant transition 
for each of the three considered S11-, -P13-, and -Di3-wave resonances. The curves are: model M2 
(full) and the meson cloud effect Eq. @ turned off (dotted) . The dashed curves correspond to the 
relevant third resonances switched off: 5n(1806) in Eq, i-i3(1893) in Ef and , and -Di3(1954) 
in E2 and ■ 

The meson cloud effect on the resonances can be illustrated by comparing the multipole 
amplitudes calculated with and without the terms within the square brackets of Eq. Q. 
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Other quantities of the coupled-channel equations are kept the same in these two calcula- 
tions. In Fig. the full curves correspond to the full M2 model, while the dotted lines are 
obtained by turning off terms within the square brackets of Eq. (jUJ), showing the importance 
of meson cloud effects in interpreting the extracted N* parameters. To further understand 
the meson cloud effects, we need to extend the present model to investigate electroproduc- 
tion data such that the Q 2 evolution of the multipole amplitude can be extracted, as has 
been done in the study of the A resonance of Ref. 0, Our effort in this direction will 
be reported elsewhere. 

In that Figure the dashed lines are obtained by switching off the relevant third reso- 
nances investigated here. These results confirm our conclusions in Sec Jlll B"3l namely, the 
Di3(1954) plays a crucial role, S , n(1806) has a significant effect, and contributions from the 
the Pi3(1893) resonance are smaller than those from the two other new resonances. 



IV. SUMMARY AND CONCLUSIONS 



The main motivation of the present work is the interpretation of recent associated 
strangeness photoproduction on the proton, which require coupled-channel formalisms. In 
the present work we have focused on the intermediate state nN, as well as the intermediate 
and final states KY interactions. 

We have first applied our formalism to the irp — > KY and KY — > KY (K = K°, K + , 
and Y = A, E°, S + ) by improving our previous work j3] and comparing successfully our 
resnfis with the existing data. We have hence fixed the factions * w and as 
well as relevant iV* parameters. Then, starting from the formalism reported in Ref. [18], we 
have developed a more advanced coupled- channel approach. For the direct jp —>■ K + A we 
have used a chiral constituent quark model j^J. The relevant data have been used to fix 
the strengths of intervening resonances within the broken SU(6) (g) 0(3) symmetry. 

Good fits to all of the available data of ir~p — > iv"°£°, 7i~p — > K°Y, , and jp — > K + A have 
been achieved. We have demonstrated that the coupled-channel effect can strongly change 
the results from the often used tree-diagram models. We have also found that the meson 
cloud effects on 'yN — > iV* are important in interpreting the extracted resonance parameters. 

This work shows that the most relevant known resonances in 7^ — > K + A process are: 
Sii(1535), Pi 3 (1900), £13(1520), and to a lesser extent Fi 5 (1680) and Fi 5 (2000). Contribu- 
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tions from three new nucleoli resonances have been extensively studied leading to convincing 
manifestations of a -D13 resonance with M = 1.954 GeV and T = 249 MeV. Rather significant 
effects due to a Su resonance with M = 1.806 GeV and T = 300 MeV is observed. A non 
negligible role is also found for a P13 resonance with M = 1.893 GeV and T = 204 MeV. 
Accounts of indications on those resonances from other sources were summarized. 

As a next step, the very new data from LEPS [f| and forthcoming polarized beam 
data from GRAAL and beam-recoil double polarization asymmetries from CLAS [69 1 ] and 
GRAAL, will hopefully clear up the experimentalsituation with respect to some inconsis- 
tencies within the present data base. Moreover, the ongoing extension of our approach to 
the 7£> — > K + Y,°, K°Y, + channels will certainly bring indeeper insights to the associated 
strangeness photoproduction processes. 

Finally, we emphasize that the present coupled-channel calculation is still far from being 
complete. While the coupling with the irN channel has been included, it is necessary to 
extend the present investigation to include the other channels, in particular the two-pion 
channels. Thus the extracted resonance parameters should be considered preliminary. But 
they could serve as the starting points for performing a more advanced coupled- channel 
calculation including additional meson-baryon channels (e.g. rjN, ujN, tittN (aN, it A, pN)) 
and to fit simultaneously all meson photoproduction data. 
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APPENDIX A: DERIVATION OF SCATTERING EQUATIONS 



In this appendix, we show that the scattering Eqs. (l)-(6) in Section III Al can be de- 
rived exactly by using the formal scattering theory given in, for example, the text book of 



Goldberger and Watson 



481 ] . We start with the following Hamiltonian 



H = Hn + V + W. 



(Al) 



where H is the free Hamiltonian, v is the non-resonant meson-baryon (MB) interaction 
with MB = jN, KY, 11N, and 



w 



rt- 



1 



(A2) 



E-H 

defines the resonant excitation by the N* — > MB vertex interaction T. 

The MB reaction amplitude T(E) is then defined 4$j by (we omit +ie in the propagator 
l/[E-H + ie}) 



or 



T(E) = (v + w)[l + 



T(E) = (v + w)[l + 
= (v + w)- 



1 



E-H 



T{E)\ 



1 



v + w)} 



E — H — v — w 

1 (E - H ) . 



E — H — v — w 
Comparing Eqs. (|A3|) and (|A5j) . we thus have 



[1 + 



-T{E)\ 



1 



E-Ho E-Ho-v-w 
We further define the non-resonant scattering matrix t by 



[E - Ho) . 



t(E) = v[l + 
= v[l + 
= [l-v 



E-Ho 
1 



t(E)} 



E-Ho-v 



1 



i-i„ 



[l+t(E) 



E — H ' 
1 



E — Hn 



Eqs. (J5ZJ and J2SJ lead to 

[1 + 



E-Ht 



■t(E)} 



E-Hq-v 
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[E - Ho] 



(A3) 

(A4) 
(A5) 

(A6) 

(A7) 
(A8) 
(A9) 
(A10) 

(All) 



and Eqs. JM1) and (lATOl) to 

[i-^r 1 l = 1+t ^E^w - (A12) 

Using Eqs. JUJ), (|A"9j) . and JHSJ), Eo.(jA"3"|) can be written as 

T{E) = [l-v-^]- 1 v + [l-v-^]- l w[l + -^—T{E)) 

III I1Q III I1Q III IIQ 

= t(E) + [1 + ((EJ-^-L—Ml + ^jfTiE)] 

= ( < B » + I 1 + WEh$ W E-H.-tv+„) < E - H "> ' (A13) 



We next use the property that 



i I .. . ( Ai4 ) 



£ - F - (v + w) ! E — H Q — v E — Hq — v 

with 



£ w = w + ^ , (A15) 

E — H n — v 



to write Eq. (|A13|) as 



T(E) = i+[1 + ((B )_^ ][ro + ro _l_ ( J_l_ [B _ Ho] 

= f + [i + 'Wfr^fe ^.^., ^ - H A ■ ( A16 > 

By using Eq. ()All|) . we then obtain 

T{E) =t+[l + t(E)-±—]t w [l + -±—t(E!)) . (A17) 

hi IIQ III IlQ 

From the separable form Eq. (|A2|) for w, it is easy to find the solution of Eq. (jA15j) 



with 



= r^-rt (A20) 
= r ^Ho rt ' (A21) 
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where 



r + = [l+t(£)-^— ^-]r+, (A22) 

f = r[l + ^±^-t{E!)) • (A23) 

Substituting Eq. ()A18jl into Eq. (IA17I) . we finally obtain 

T{E) = t + r + ^_^_ s r . (A24) 

Taking the matrix elements of the relevant equations given above between the channels 
a,b,c = '-yNjirN, KY and noting that Hq\N* >= M^t*\N* > in the center of mass frame, 

i 

G a =< a\ e ^ Hq [a >, and T N * <a =< N*\T\a >, we then obtain Eqs. (l)-(6) in Section Hi Al 
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